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Abstract: The Riodinidae is one of the lepidopteran butterfly families. This study describes the complete mitochondrial genome of 
the butterfly species Abisara fylloides, the first mitochondrial genome of the Riodinidae family. The results show that the entire 
mitochondrial genome of A. fylloides is 15301 bp in length, and contains 13 protein-coding genes, 2 ribosomal RNA genes, 22 
transfer RNA genes and a 423 bp A+T-rich region. The gene content, orientation and order are identical to the majority of other 
lepidopteran insects. Phylogenetic reconstruction was conducted using the concatenated 13 protein-coding gene (PCG) sequences of 
19 available butterfly species covering all the five butterfly families (Papilionidae, Nymphalidae, Peridae, Lycaenidae and 
Riodinidae). Both maximum likelihood and Bayesian inference analyses highly supported the monophyly of Lycaenidae+Riodinidae, 
which was standing as the sister of Nymphalidae. In addition, we propose that the riodinids be categorized into the family Lycaenidae 


as a subfamilial taxon. 
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The typical metazoan mitochondrial genome 
(mitogenome) contains 37 genes, including 13 protein- 
coding genes (PCGs), 2 rRNA genes and 22 tRNA genes, 
and a non-coding area (i.e., the control region or the 
A+T-rich region) (Wolstenholme, 1992; Boore, 1999). 
Maternal inheritance, lack of recombination and an 
accelerated evolutionary rate compared with the nuclear 
genome have all contributed to the increased use of 
mitogenomes, which is one of the key methods in fields 
such as phylogenetics, comparative and evolutionary 
genomics, molecular evolution and population genetics 
(Ballard & Whitlock, 2004; Simonsen et al, 2006). At 
present, mitochondrial genomes have already been 
determined in a variety of insect groups covering nearly 
200 species. However, reported complete mitogenomes 
are relatively scarce for lepidopterans and especially for 
butterflies. To our knowledge, as of October 2012 only 
about 20 butterfly species covering five butterfly families 
(Table 2) have been reported or deposited into the 
GenBank, but only one butterfly family, the Riodinidae, 
still lack corresponding data. 

The phylogenetic position and taxonomic ranking of 


Science Press 


the butterfly family Riodinidae among butterfly lineages 
are still controversial issues among entomologists. Some 
scholars suggest that the riodinids are closely related to the 
lycaenids, considering the similarities in morpholo gical 
character, behavior, and host plants between the two (slug- 
like larvae, pupa contigua, ants associated) (Ackery, 1984; 
Chou, 1998; de Jong et al, 1996; Ehrlich, 1958; Scott, 
1985). Moreover, they are usually classified into the 
Lycaenidae family as a subfamilial taxon. Some consider 
the riodinids a unique family parallel to the Lycaenidae 
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family (Campbell et al, 2000; de Jong et al, 1996; 
Kristensen, 1976; Shou et al, 2006), and others propose 
that the riodinids are more closely related to the 
nymphalids than to the other butterfly groups, such as the 
lycaenids (Martin & Pashley, 1992; Robbins, 1987, 1988). 
This study sequenced the first complete 
mitochondrial genome of the Abisara fylloides, a 
representative species of the family Riodinidae, by long 
PCR and primer walking techniques. Its genetic structure 
was preliminarily compared with those of other available 
butterfly species. The maximum likelihood (ML) and 
Bayesian inference (BI) phylogenetic trees of the A. 
fylloides and other available butterfly representative 
species were reconstructed based on concatenated DNA 
sequences of the 13 protein-coding genes (PCGs), and aim 
to clarify their phylogenetic relationships and further 
provide new information about the structure, organization 
and molecular evolution of the lepidopteran mitogenomes. 


MATERIALS AND METHODS 


Sample collection and DNA extraction 
Adult individuals of A. fylloides were collected in 


Jinghong, Yunnan Province, China in August 2006. After 
collection, sample specimens were preserved in 100% 
ethanol immediately and stored at —20 °C before DNA 
extraction (specimen No. ZHF07). Total Genomic DNA 
was isolated using the proteinase K-SiO, method as 
described by Hao et al (2005). 
PCR amplification and sequence determining 

The multiple sequence alignments were conducted 
using the software Clustal X 1.8 based on the mitogenome 
sequences of Coreana raphaelis, Artogeia melete, Troides 
aeacus available from GenBank and those of Argyreus 
hyperbius, Acraea issoria, Calinaga davidis, Pieris rapae 
determined in our laboratory (Thompson et al, 1997). The 
long PCR primers, which may cover the whole 
mitogenome, were designed according to the conserved 
regions by the software Primer premier 5.0 (Singh et al, 
1998) (Table 1). Seven short fragment sequences (500-700 
bp) of сохі, cox2, cox3, cytb, пайї, rrnL and rmS were 
amplified using insect universal primers (Caterino & 
Sperling, 1999; Simmons & Weller, 2001; Simon et al, 
1994). All the primers were synthesized by the Shanghai 
Sheng gong Biotechnology Co. Ltd. 


Table 1 List of PCR primers used in this study 





Genes Forward primers (5) Reverse primers (3’) ы C) 
сох1* GGTCAACAAATCATAAAGATATTG TAAACTTCAGGGTGACCAAAAAT 50.0 
сох1-сох2 ТТАТТТСТАТСАССССТАС ATAGCAGG AAGATTGTTC 47.5 
сох2* GAGACCATTACTTGCTTTCAGTCACT CTAATATGGCAGATTATATGTATGG 49.5 
cox2-cox3 TTTTATTGCTCTTCCATCT TTATTCCTCATCGTAATCC 48.5 
cox3* TATTTCAATGATGACGAGAT CAAATCCAAAATGGTGAGT 49.8 
cox3-nad5 TTTATAGCAACAGGATTTC CATCAACTGGTTTAACTTT 45.5 
nad5 AAAACTTCCAGAAAATAATCTC TTGCTTTATCTACTTTAAGACA 46.5 
nad5-cytb AATTATACCAGCACATAT TTATCGACTGCAAATC 47.1 
cytb* TATGTACTACCATGAGGACAAATAT ATTACACCTCCTAATTTATTAGGAAT 47.0 
cytb-nadl TCCTGCTAACCCTTTAGTCA AGGTAGATTACGGGCTGTT 48.0 
nad1* CGTAAAGTCCTAGGTTATATTCAGATCG ATCAAAAGGAGCTCGATTAGTTTC 52.0 
nadl-rrnL AGCCCGTAATCTACCTAA TAAGACGAGAAGACCCTAT 47.0 
rmL* CGCCTGTTATCAAAAACAT CCGGTCTGAACTCAGAT 45.5 
rrnL-rrnS AGACTATTGATTATGCTACCT TAAGAATCTAATGGATTACAA 46.5 
rm s* CTTCTACTTTGTTACGACTTA T AATTTTGTGCCAGCAGTTG 50.0 
rrnS-nad2 AGAGGGTATCTAATCCGAGTTT TGGCTGAGAATTAAGCGATA 49.5 
nad2-coxl ATACAGAAGCAGCATTA AGAAGGAGGAAGAAGTCAA 52.0 


*: universal primer. 


Seven partial gene sequences were initially 
sequenced under the following conditions: an initial 
denaturation at 94 °С for 5 minutes, then denaturation at 
94 °С for 1 minute for a total of 35 cycles; annealing at 
45—55 °С for 1 minute and extension at 72 °С for 2 
minutes plus 30 seconds; final extension at 72 °С for 10 
minutes. Long PCRs were performed using TaKaRa LA 
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Taq polymerase with the following cycling parameters: 
an initial denaturation for 5 minutes at 95 ?C; followed 
by 30 cycles at 95 °С for 55 seconds, 45—55 °С for 2 
minutes, 68 ?C for 2 min and 30 seconds; and a 
subsequent final extension step of 68 °С for 10 minutes. 
The PCR products were separated by electrophoresis 
in a 1.296 agarose gel and purified using the DNA gel 
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extraction kit (TaKaRa). All PCR fragments were 
sequenced directly after purification with the QIA quick 
PCR Purification Kit reagents (QIAGEN). Internal primers 
were applied to complete sequences by primer walking 
(detailed primer information will be provided upon request). 
АП fragments were sequenced for both strands. 


Data analysis 

We used DNASIS MAX (Hitachi) for sequence 
assembly and annotation. Protein-coding genes and 
rRNA genes were identified by sequence comparison 
with other available insect mitochondrial sequences. The 
tRNAs were identified by tRNAscan-SE v.1.21 (Lowe & 
Eddy, 1997). The putative tRNAs, which were not found 
by tRNAscan-SE, were identified by a sequence 
comparison of A. fylloides with the other lepidopteran 
tRNAs. PCGs were aligned with the other available 
lepidopteran mitogenomes using DAMBE software (Xia 
& Xie, 2001). The tandem repeats in the A+T-rich region 
were predicted using the Tandem Repeats Finder online 
(http://tandem.bu.edu/trf/trf.html) ^ (Benson, 1999). 
Nucleotide composition was calculated using PAUP 
4.0b10 (Swofford, 2002). The mitogenome sequence 
data have been deposited in GenBank under the 
accession number HQ259069. 


Phylogenetic analysis 

Phylogenetic analyses were performed on 19 
representative species including A. fylloides, covering all 
the six families of butterflies. The multiple aligning of 
the concatenated nucleotide sequences of the 13 
mitochondrial PCGs of the 19 species (Table 2) was 
conducted using ClustalX 1.8. The phylogenetic trees 
were reconstructed with the maximum likelihood (ML) 
and Bayesian inference (BI) methods, using the moth 
species Adoxophyes honmai (GenBank accession 
number of mitogenome: NC014295) as the outgroup. In 
both phylogenetic analyses, the third codon position of 
all the sequences was excluded. The ML analyses were 
conducted in PAUP 4.0b10 by using TBR branch 
swapping (10 random addition sequences) as a search 
method. The model GTR+I+T was selected as the best fit 
model using Modeltest 3.06 (Posada & Crandall, 1998) 
under the AIC scores, and the bootstrap values of the ML 
tree were evaluated via the bootstrap test with 1 000 
iterations. The Bayesian analysis was performed using 
MrBayes 3.1.2 (Huelsenbeck & Ronquist, 2001) with the 
partitioned strategy (13 partitions: cox], cox2, cox3, atp8, 
atp6, nad1, nad2, nad4, nad4L, nad5, nad6 and cytb), and 
the best substitution model for each partition was 
selected as in the ML analysis. The MCMC analyses 
(with random starting trees) were run with one cold and 
three heated chains simultaneously for 1,000,000 
generations sampled every 100 generations with a burn- 
in of 25% until the average standard deviation of split 
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frequencies to be less than 0.01, which means that 
convergence was reached. 


RESULTS AND DISCUSSION 


Genome structure and organization 

The complete mitogenome is 15 301 bp in length, 
encodes 37 genes in all. It contains 13 protein, 22 tRNA, 
2 rRNA genes and a non-coding high A+T content region 
(Figure 1, Table 3). Its structure and organization are 
identical to those of the majority of other lepidopterans 
(Bae et al, 2004; Cha et al, 2007; Cameron & Whiting, 
2008; Hao et al, 2012; Hong et al, 2008, 2009; Hu et al, 
2010; Kim et al, 2009b; Ji et al, 2012; Junqueira et al, 
2004; Wang et al, 2011), though a few lepidopterans, 
such as three Thitarodes species, were reported to 
possess the ancestral gene order trnI-trnQ-trnM instead 
of the trnM-trnI-trnQ (Cao et al, 2012) . 

Eight overlapping sequences totaling 61 bp are 
located throughout the A. fylloides mitogenome, with 
size ranging from 2 to 35 bp, of which the longest (35 bp) 
is located between the cox2 and the tRNA genes. In 
addition, 17 intergenic spacers ranging from 1 to 45 bp in 
length are found in the mitogenome. Among these 
spacers, the longest is located between the КМА" and 
nad2 genes, the other 16 spacers are scattered throughout 
the whole genome (Figure 1, Table 3). Most of these 
spacer regions are arranged relatively compactly 


Abisara fyuoides 


mtDNA 15301 bp 





Figure 1 Circular map of the mitochondrial genome of 
Abisara fylloides 

Gene names not underlined indicate the direction of transcription from left to 

right and those underlined indicate right to left. Transfer RNA genes encoded 

by H and L strands are shown outside and inside the circular gene map, 

respectively. Transfer RNA genes are indicated by the IUPAC-IUB single letter 

amino acid codes, while L1, L2, S1, S2 represent tRNA-Leu(UUR), tRNA- 


Leu(CUN), tRNA-Ser(AGN) and tRNA-Ser(UCN), respectively. 
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Table3 Organization of the Abisara fylloides mitochondrial genome 





Gene Direction Position Size (bp) Intergenic length" Start codon Stop codon 
tRNAM* F 1-67 67 0 
tRNA! F 68-131 64 -3 
tRNA?" R 129-197 69 45 
nad2 F 243-1 256 1014 1 ATT TAA 
tRNA"? F 1 258-1 324 67 -8 
tRNA R 1317-1382 66 3 
tRNA’ R 1386-1 451 66 6 
сох1 Е 1 458-2 989 1 532 3 CGA T 
ВМА 0900 Е 2 993-3 059 67 0 
cox2 F 3 060-3 738 679 -3 ATG T 
tRNA F 3 736-3 806 71 
tRNA‘? F 3 807-3 871 65 
atp8 F 3 872-4 033 162 -7 ATC TAA 
atp6 F 4 027-4 714 688 -2 ATG TAA 
сох3 Е 4 713-5 501 789 2 ATG TAA 
tRNAGY F 5 504—5 569 66 0 
nad3 F 5 570—5 923 354 3 ATT TAA 
tRNAA" F 5 927—5 997 71 3 
tRNA^'? F 6 001-6 064 64 1 
tRNA n F 6 066—6 131 66 16 
ЕТА АЧА» Е 6 148—6 208 61 1 
tRNAS" F 6 210-6 275 66 -2 
tRNA’ R 6 274-6 340 67 0 
nad5 R 6 341-8 108 1 768 ATT T 
tRNA" R 8 109-8 179 71 
nad4 R 8186-9 518 1333 E ATG T 
nad4L R 9 517—9 805 289 2 ATG TAA 
tRNA!” F 9 808—9 870 63 0 
tRNA"? R 9 871-9 937 67 5 
nad6 F 9 943-10 467 525 3 ATA TAA 
cytb F 10 471-11 622 1152 -2 ATG TAA 
tRNASEHUCN) F 11 621-11 685 65 17 
пайї R 11 703-12 641 939 1 ATG TAA 
tRNAP CUN R 12 643-12 710 68 0 
rrnL R 12 711-14 044 1334 0 
tRNA I R 14 045-14 107 63 0 
rrnS R 14 108-14 878 771 0 
A+T-rich region 14 879-15 301 423 


`: In the column intergenic length, the positive number indicates interval base pairs between genes, while the negative number indicates the overlapping base 


pairs between genes. 


compared with other insect mitogenomes (Cameron & 
Whiting, 2008; Hao et al, 2012; Hong et al, 2009; Hu et 
al, 2010; Kim et al, 2009b; Ji et al, 2012; Junqueira et al, 
2004; Wang et al, 2011). 

The A. fylloides A+T-rich region is flanked on one 
side by the rmS and on the other side by the tRNA Y° 
genes. This region exhibits the highest A+T content 
(91.096) (Table 4), and spans 423 bp (Table 3, Table 4). A 
sequence analysis of the A+T-rich region revealed that it 
contained some structures typical of other lepidopteran 
mitogenomes: (1) at 20 bp downstream of the small 
subunit rRNA gene, there is a structure including a 
motif ‘АТАСА? which is very well conserved in all 
sequenced lepidopteran insects, and a 18-bp polyT 
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stretch, both of which have been suggested as the origin 
of minority or light strand replication (On) and to play a 
regulatory role (Kim et al, 2009a; Lutz-Bonengel et al, 
2004; Saitou et al, 2005; Yukuhiro et al, 2002). (2) 
Between the sites 15 151 and 15 168 there is a 
microsatellite-like (TA) element, which is also found 
in the majority of other lepidopterans; (3) There is 
another motif “ATTTA” of unknown function located 
from 15 139 to 15 143 upstream of ће (TA)s, which is 
also typical of the other lepidopterans. In addition, to 
our great surprise, an unexpected short microsatellite- 
like repeating region (ТА), downstream of the (ТА) 
was detected in the А+Т-гїсһ region, and this has not 
been reported in any other lepidopterans. 
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Table 4 Nucleotide composition and skewness in different regions of the Abisara fylloides mitogenome 


Nucleotide composition (%) 








Region в AT-skew GC-skew 
(bp) T б А G AKT 

Whole genome 15301 41.7 11.3 39.5 7.5 812 —0.027 —0.202 
Major-strand PCGs 6 927 44.5 12.1 34.1 9.3 78.6 —0.132 —0.131 
Minor-strand PCGs 4329 47.6 6.3 34.1 12.0 81.7 —0.165 0.311 
Whole PCGs 11 224 45.8 9.9 34.0 10.3 79.8 —0.148 0.020 
Whole tRNA 1 461 40.0 TA 41.5 10.7 81.5 0.018 0.163 
rmL 1 334 38.4 48 47.0 9.8 85.4 0.101 0.342 

rms 771 39.9 4.9 45.7 9.5 85.6 0.068 0.319 
A+T-rich region 423 48.9 4.5 42.1 4.5 91.0 —0.075 0.000 


Base composition bias 

The base composition of the А. fylloides 
mitogenome shows an A and T bias (Table 4). The 
whole А+Т content of the mitogenome 15 up to 81.295, 
ranging from that of Argyreus hyperbius (80.81%) to 
that of Coreana raphaelis (82.66%). Like most other 
metazoan mitogenomes, the A+T content of the A+T- 
rich region, which is located between the rrnS and 
tRNA М“ genes, is the highest (91.0%) in all known 
butterfly species except for P. rapae (Pieridae) to date. 
The base contents of A, T, C, G are 39.5%, 41.7%, 
11.3%, 7.5%, respectively, indicating a relatively 
higher A+T content (81.2%). These phenomena 
commonly exist in the protein-coding genes, which 
have a relatively lower A+T content (79.8%), and the 
tRNA and rRNA genes in insects. 


Protein-coding genes 

The sequences of the 13 Abisara fylloides PCGs 
are 11 224 bp in length, including 3 730 codons 
(excluding termination codons). Twelve of the 13 
PCGs use standard ATN as their start codon except for 
the coxl gene, and eight of these 12 PCGs begin with 
ATA or ATG (Methionine) and the other four begin 
with АТТ or ATC (Isoleucine) (Table 3). The start 
codons for coxl gene of lepidopteran insects have 
usually been a controversial issue. In general, there is 
no typical start codon, especially for the сох1 gene, 
which usually uses CAG (R) as start codon in insects, 
including the lepidopterans. However, some scholars 
reported unusual start codons, such as the 
trinucleotide TTG (Bae et al, 2004; Hong et al, 2008), 
ACG (Lutz-Bonengel et al, 2004), GCG (Nardi et al, 
2003), the tetranucleotide ATAA, ATCA and ATTA 
(Clary & Wolstenholme, 1983; de Bruijn,1983; Kim et 
al, 2006), and the hexanucleotides TATTAG (Flook et 
al, 1995), TTTTAG (Yukuhiro et al, 2002), TATCTA 
(Coates et al, 2005), ATTTAA (Beard et al, 1993; 
Mitchell et al, 1993) for coxl in some other insect 
species. However, the CGA is present as a conserved 
region for all lepidopteran insects reported, and thus 
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we tend to consider that CGA is the сох1 start codon 
for A. fylloides as Kim et al (2009b) suggested. As for 
stop codons, 9 of the 13 PCGs use standard TAA 
except for the coxl, cox2, nad4 and nad5 genes, all of 
which terminate at a single thymine (Table 3), and this 
case was also found in all other lepidopterans reported 
to date. For more details on this phenomenon, please 
refer to the discussions of Kim et al. (2010). 

The PCG amino acid sequence variation analysis 
showed that there were 3 755 homologous sites, of 
which 1 964 are conserved, 1 791 are variable, and 1 
128 are parsimony informative. Among the twenty 
amino acids in the 13 PCGs, six (Leu, Met, Ile, Phe, 
Asn, and Tyr) were used more frequently than the 
others, and their usage frequencies were higher than 
the average, whereas the other 14 amino acids less 
used (Table 5). 


rRNA and tRNA genes 

The large subunit rRNA and small subunit rRNA 
genes of the A. fylloides are 1 334 and 771 bp in length, 
respectively. As in other lepidopterans, these two genes 
are located between tRNA!” and tRNA, and 
between tRNA I and A+T-rich region respectively 
(Cameron & Whiting, 2008; Hao et al, 2012; Hu et al, 
2010; Kim et al, 2009a; Salvato et al, 2008; Wang et al, 
2011; Yang et al, 2009). 

The A. fylloides mitogenome harbors 22 tRNA 
genes, which are scattered throughout the whole genome 
and ranged in length from 61 to 71 bp. Except for the 
tRNA SGN, which lacks the DHU loop, all tRNAs are 
shown to be folded into the cloverleaf secondary 
structures, within which all amino acid acceptor stems 
have 7 base pairs, and all anticodon stems have 5 base 
pairs. This was also found in all the other lepidopterans 
determined to date (Cameron & Whiting, 2008; Hao et al, 
2012; Hu et al, 2010; Kim et al, 2009a). 

A total of 27 pairs of base mismatches were 
detected in all the predicted tRNA secondary structures, 
among which 17 are GU, 6 are UU, 2 are AA, and 2 are 
AC. The AA mismatches occur at the anticodon stem of 
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Table 5 Codon usage of the protein-coding genes of the Abisara fylloides mitogenomes 








Codon N (RSCU) Codon N (RSCU) Codon N (RSCU) Codon N (RSCU) 
UUU (F) 356 (1.86) UCU (S) 108 (2.59) UAU (Y) 175 (1.83) UGU (С) 28 (1.81) 
UUC (F) 26 (0.14) UCC (S) 16 (0.38) UAC (Y) 16 (0.17) UGC (C) 3 (0.19) 
UUA (L) 488 (5.30) UCA (S) 94 (2.25) UAA (*) 0 (0) UGA (W) 84 (1.87) 
UUG (L) 12 (0.13) UCG (S) 2 (0.05) UAG (*) 0 (0) UGG (W) 6 (0.13) 
CUU (L) 34 (0.37) CCU (P) 70 (2.33) CAU (H) 62 (1.85) CGU (R) 14 (1.04) 
CUC (L) 3 (0.03) CCC (P) 9 (0.3) CAC (H) 5 (0.15) СОС (В) 0 (0) 
CUA (L) 14 (0.15) CCA (P) 38 (1.27) CAA (Q) 57 (1.93) CGA (R) 37 (2.74) 
CUG (L) 1 (0.01) CCG (P) 3 (0.1) CAG (Q) 2 (0.07) CGG (R) 3 (0.22) 
AUU (D 424 (1.88) ACU (T) 74 (2.1) AAU (N) 251 (1.81) AGU (S) 25 (0.6) 
AUC (D 27 (0.12) ACC (T) 7 (0.2) AAC (N) 26 (0.19) AGC (S) 0 (0) 
AUA (M) 283 (1.81) ACA (T) 59 (1.67) AAA (К) 97 (1.92) AGA (S) 89 (2.13) 
AUG (M) 29 (0.19) ACG (T) 1 (0.03) AAG (K) 4 (0.08) AGG (S) 0 (0) 
GUU (V) 74 (2.39) GCU (A) 73 (2.52) GAU (D) 55 (1.83) GGU (G) 56 (1.17) 
GUC (V) 4 (0.13) GCC (A) 8 (0.28) GAC (D) 5 (0.17) GGC (G) 1 (0.02) 
GUA (V) 45 (1.45) GCA (A) 33 (1.14) GAA (E) 65 (1.78) GGA (G) 100 (2.08) 
GUG (V) 1 (0.03) GCG (A) 2 (0.07) GAG (E) 8 (0.22) GGG (G) 35 (0.73) 





n: frequency of codon used; RSCU: relative synonymous codon usage; *:stop codon. 


Stop codons were excluded in total codon counts. 


tRNA'^ and the TyC loop of tRNA"; the AC 
mismatches occur at the amino acid acceptor stem of 
tRNA " and the anti codon stem of tRNA"), Among 
the 6 UU mismatches, two occur at the amino acid 
acceptor stems of tRNA г and tRNA F, two at the 
anti codon stems of {ВМА (ЧУ and ВМА", and 
another two at the anti codon stem of tRNA?* CV, 
respectively (Figure 2). 


Phylogenetic analysis 

At present, for the phylogenetic positions of the 
riodinids within —papilionid butterflies, most 
morphological studies place them as most closely related 
to the lycaenids and identify the nymphalids as the 
closest relatives to this riodinid+lycaenid clade (de Jong 
et al, 1996; Ehrlich & Ehrlich, 1967; Kristensen, 1976; 
Scott & Wright, 1990). These relationships have been 
inferred using a variety of phylogenetic methods and are 
supported by a number of adult, larval and pupal 
synapomorphies. ^ Additionally, molecular (DNA 
sequence of the mitochondrial МАРНІ gene) or 
molecular plus morphological evidence also result in a 
monophyletic interpretation of the Riodinidaet+ 
Lycaenidae, and their sister relationship to the 
Nymphalidae (Wahlberg et al, 2005; Weller et al, 1996). 
However, based on a cladistic analysis of four foreleg 
characters with nine character states, Robbins (1988) 
suggested that the Riodinidae are more closely related to 
the Nymphalidae than to the Lycaenidae, and this result 
is supported by the nuclear 28S rRNA gene sequence 
data (Martin & Pashley, 1992). 

There are two opinions regarding the taxonomic 
rank of riodinids. First, some previous studies suggest 
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that the riodinids should be categorized into the family 
lycaenids as a subfamilial taxon in light of their 
morphological characters (Chou, 1998; de Jong et al, 
1996; Ehrlich, 1958; Kristensen, 1976; Scott & Wright, 
1990), and this opinion is supported by the molecular 
studies of Zou et al (2009) and Hao et al (2007). Other 
studies postulate that the riodinids should be classified as 
a separate family parallel to Lycaenidae (Harvey, 1987; 
Martin & Pashley, 1992; Robbins, 1988; Weller et al, 
1996), and this view is supported by molecular 
phylogenetic studies based on data from the wingless 
gene by Campbell et al (2000) and the combined analysis 
of sequences of the nuclear Ef-la, wingless, and 
mitochondrial COI genes by Wahlberg et al (2005). 

The ML and Bayesian trees of this study (Figure 3) 
showed that all the butterfly taxa in this study did not 
form a monophyletic unit, and a similar case was 
reported by Hao et al (2012). Nonetheless, both the ML 
and BI trees indicated that all the butterfly species were 
grouped into five distinct lineages: 1) the Papilionidae, 
including papilionids and parnassids; 2) Hesperiidae; 3) 
Pieridae; 4) Nymphalidae; 5) Lycaenidae+Riodinidae. 
The monophyly of Lycanidae + Riodinidae was strongly 
supported with a 100% bootstrap value in ML, and with 
1.00 posterior probability value in BI. Thus, considering 
the results of Heikkilà et al (2012), which indicate the 
monophylies of lycaenids and riodinids, it is reasonable 
to propose that the two groups may be sisters, though the 
taxa sampling of riodinids in this analysis is extremely 
limited. Additionally, based their congruent genetic 
divergences compared with those between other butterfly 
subfamilies, the riodinids should be categorized into the 
Lycaenidae family as a subfamilial taxon. 
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Figure2  Predicated clover-leaf secondary structures for the mitochondrial tRNA genes of Abisara fylloides 
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Figure3 The Bayesian inference (BI) and maximum likelihood (ML) phylogenetic trees of main butterfly lineages based on 13 


protein-coding gene sequences (Numbers on each node correspond to the posterior probability values of the BI analysis and the ML 


bootstrap percentage values for 1 000 replicates of ML analysis) 
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